Does negative refraction make a perfect lens? 
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Abstract 

A discussion of a question, studied earlier by V.Veselago in 1967 and by J. 
Pendry in 2000, is given. The question is: can a slab of the material with 
negative refraction make a perfect lens? Pendry's conclusion was: yes, it can. 
Our conclusion is: no, in practice it cannot, because of the fluctuations of 
the refraction coefficient of the slab. Resolution ability of linear isoplanatic 
■ optical instruments is discussed. 
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£5 : 1 Introduction 

p 

Negative refraction has been studied extensively (see pQ and references therein). In 
[I] the following question is discussed. Assume that a plane electromagnetic (EM) 
wave E = e2e t ( k ^- x + k 3 z ) := e 2 e tklX u (z), u (z) = e lk3Z , is incident from the region 
z < onto an infinite slab (layer) < z < d, filled with a material with negative 
refraction (negative index). This means that e < and \x < in the slab, where e 
^ ■ and fx are constant dielectric and magnetic parameters. The refractive index 



n = (^ = (| e ||^ = -(| E || /t r<0. 

The time dependence is given by the factor e~ tuJt and this factor is omitted through- 
out. By co > the frequency is denoted. The vectors {ej}? =1 form a Cartesian 
orthonormal basis of M 3 . Outside the slab, in the region z < and z > d, we as- 
sume that e = eo and (i = /jlq, where 6q — fJ^o — 1. It is assumed in [J| that in the slab 
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e = — 1 and fi = — 1, so that n = — 1 in the region < z < d. By fcj the components 
of the wave vector k are denoted, k = k\e\ + k^e^, \k\ 2 = \ki\ 2 + \k 3 \ 2 = u 2 efi. In our 
case &2 = 0, and e/i > 0. Let us denote the constants e and \x in the slab by t\ and 
Mi- 

The governing equations are the Maxwell's equations: 

V x E = iupiH, V x H = -icueE, (1) 

and the tangential components of the electric and magnetic fields E t and H t are 
continuous across the interfaces z = and z = d. 
Our calculations prove the following conclusion: 

If n = —1 in the slab, then the image, of a point source, located at the point 
(0,0, -/), is a point (0,0, 2d - /). 

This conclusion agrees with the one obtained earlier in [13] and [1]. In [1] only 
the propagation of plane waves is discussed. However, the point source is a linear 
combination of the plane waves (including the evanescent waves), so there is no loss 
of generality in considering propagation of the plane waves. For example, a well 
known representation of a point source as a linear combination of the plane waves 
is the formula: 



e ik\x\ ^ f f 00 Qi{kix\+k2X2+x 3 ^/k 2 —k\—k 2 l ) 

- — |— r = / / dkidko, 

4k\x\ 2niJJ_ 00 ^k 2 - k 2 - k 2 

where X\ = x, X2 = y, £3 = z, and \/k 2 — k\ — k\ > if k\+k\ < k 2 , \Jk 2 — k\ — k\ = 
i^\k 2 -k 2 -k 2 \ if k\ + k\ > k 2 . 

Let us derive the formulas for the electromagnetic field, propagating through the 
slab, and prove the conclusion, formulated above. 

Let us look for a solution to flTJ of the form: 

e ik 3 z _|_ re -ik 3 z^ Z <0, 

E = e 2 e iklX u(z), u{z) = { Ae ik * z + Be- ik * z , < z < d, (2) 

te ik3Z , z> d, 

where the four coefficients r, A, B, t, are to be found from the four boundary condi- 
tions: 



Here 



1 +r = A + 


B, h(l-r) 


¥ 


(3) 


Ap' + - = 
pi 


tp, ^ [Ap' 
Mi V 


B \ 

- — = thp. 
p'J 


(4) 


1 

V 


:= e ifc 3 d , p := 


ik 3 d 


(5) 
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k , . = i v / io 2 e 1 fi 1 -kf>0, k 2 < Lu 2 ein u . . 

3 ' \ i^\u 2 e 1 n 1 - kj~\, kl>uj 2 e l n u 

and the boundary conditions ([3]) and (j3J) come from the continuity of the tangential 
components of the electric and magnetic fields, that is, [e%, E] and [e^, H] = ^ 3 ^^ ; 
across the planes z = and z = d, respectively. Here [e, E] denotes the vector 
product of two vectors. The four equations Q for the four unknowns r, A, B, t, 
can be solved analytically: 

r = 7 l ~P' 2 t = P -- l2 ~ l (7) 

'9/9 1 ' 9 /9 1 ' \ / 



1 + b ~f 2 p' 2 - 1 ' 2 7V 2 - 1 ' 

where p' and p are defined in dSJ) , and 

7-=— 6-=^ (9) 
7 ' 1 + 6' M 3 ' [) 

If ei = yUi = —1, then k' 3 = k 3 , b = —1, 7 = 00, 

r = 0, t = e" 2ifc3d , A = 0, 5 = 1. (10) 

The coefficients r and t are called the reflection and transmission coefficients, re- 
spectively. In formulas (12), (21) in [4J by T and T s not the standard transmission 
coefficients are denoted, but the values u(z) at z — d. The values of the field 
everywhere in the space are given by formulas (j2J), ©-©• 



2 A discussion of the exact solution 

Let us draw some physical conclusions from the the basic formulas ([2]), ((SJ)-(jSD- 
First, consider the simplest case e\ = \i\ — — 1. Then formulas (j2J), (|T0l) hold. These 
formulas imply that a ray, passing through a point (0, 0, — /), / > 0, with parameters 
k\ > 0, k 3 > 0, hi = 0, will have parameters (ki, 0,-/03) after passing through 
the interface z = 0, and again the parameters (A;i,0, ks) after passing through the 
interface z = d. In this geometrical argument it is assumed that d > f. The ray is 
the straight line orthogonal to the wave front of the plane wave. The ray issued from 
the point (0,0,—/) with parameters (ki,0,k 3 ), hits the plane z = at the point 
(/to, 0, 0), where h = /tan# > tan6> = jg, and 6*o is the acute angle that the ray 
makes with z-axis. The z-axis is perpendicular to the slab. After passing the plane 
z = 0, the ray hits the z-axis at the point (0, 0, /) and then it hits the line z = d, 
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y = 0, at the point (—hi, 0, 0), where hi = (d — f) tan# - Finally, after passing the 
plane z = d, the ray hits the z-axis at the point (0, 0, 2d — f). 
The conclusion is : 

Regardless of the size of 8 , all the rays, passing through the point (0,0,—/), 
intersect the z-axis at the point 2d — f , d > f > 0. 

This means that the slab with ex = /ii = — 1 acts like a "perfect lens", as is 
claimed in [13] and [I] . There is no need to discuss separately the case of evanescent 
waves, i.e., the case when uj 2 eifii < h\. Our derivations are different from those in 
[4], but the conclusion in the case ei — fii = — 1 is the same. 

However, in [I] he following questions are not discussed: 

1) One cannot have e = — 1 and /i = — 1 exactly due to defects in the material 
of the slab. Will the above conclusion, (namely, that the image of a bright point 
F := (0, 0, — /) is a point (0, 0, 2d — f), < / < d), remain true if k' 3 ^ k 3 ? 

2) Suppose that the slab is not infinite. Will the diffraction from the boundary of 
the slab invalidate the principal conclusion? 

We will not discuss the second question in detail, and restrict ourselves to the 
following argument. If the size of the slab, which we defined as a maximal radius R 
of a cylinder, inscribed into the slab, is much larger than the wavelength Ao = 
ko = oj-s/eofiQ, R ^> A and R^$> I, where I is the size of the object, then one may 
neglect the diffraction of light at the boundary of the slab, provided that the object 
is located near z-axis, far away from the boundary of the slab. 

The first question we discuss in more detail. The aim of our discussion is to 
conclude that if e\ or \ii differ slightly from —1, so that n is not equal to —1 exactly, 
then the image of a point (0, 0, — /) is no longer a point. More precisely, not all the 
rays, passing throught the point F := (0,0, — /), pass through the common point of 
the z-axis after the plane z = d. It turns out that the point of itersection of the ray, 
passing through F at the angle 9 with the z-axis in the region z > d, will depend on 
6\) if \k' 3 \ ^ \k 3 \. Assume, for example, that fa — — 1 + 8^, ei — — 1 + 5 e , 5^ = 8 e = 5, 
where 5 > is a small number. Then 

K = ^(-l + cyc-i + cy-fci 2 « sJ^-ki-^(5, + 5 e ) 




Here we have neglected the term of order 0(S 2 ), and assumed that \k 3 \ ^> 5 > 0. 
Thus § = 1 - 8i, where 5i := ^. 

If k?j = 0(5) then j^l = 0(5 1 ^ 2 ). In this case || can take rather arbitrary values 
depending on the ratio r^CT , and therefore the deviation of the corresponding rays 
from the point (0, 0, 2d — f) may be much greater than 0(5). 
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If |&3 1 3> S, then we repeat the geometrical arguments given for the case k' 3 = k^, 
and get the following conclusions: the ray, passing through the point F := (0, 0, — /) 
and having angle 9q with the z-axis, hits the plane z = at the point (ho, 0,0), 
h = /tan# ? then it hits the z-axis at the point (0,0,/'), where tan 9' = y- = j^, 
so /' = h cot9', and then it hits the plane z = d at the point (—^,0,0), where 

h' h 

-j-^jj = tan# = y- Finally, after passing the plane z = d, this ray hits the z-axis at 
the point (0, 0,d + z'), where z' = h\ cot 6q. Therefore, 

d + z = d + (d - f tan O cot 9') tan 9' cot 9 = d(l + tan 9' cot 0„) - /. 

One has tan6> = f^, tan6>' = |f, so tan 0' cot #0 = |f = j^-. 

Conclusion: T/ie ray, passing through the point F = (0,0,— f) under the angle 
9 with z-axis, hits the z-axis, after passing the slab, at the point whose coordinate 
(0, 0, d — / + dtan^' cot 9q) depends on the angle 9q. Therefore, the image of a point 
source, located at the point F, will not be a point. 

At a fixed uj > the deviation of the image of a bright point depends on the 
deviation of |f from 1. This implies a possibility to make changes in e\, if it is prac- 
tically easier than to make changes in fii, in order to compensate for the deviations 
of /ii from —1. 

3 On the resolution ability of linear isoplanatic 
optical instruments 

Following L. Mandel'shtam [2], assume that a linear isoplanatic optical instrument 
S creates an image of an object as follows. A bright point (a point source), is 
located on a plane Pi at a point y G M 3 , and the focal plane is located to the left 
of S. The bright point produces an image h(x, y) = h(x — y) on the image plane 
P%, x G Pi, and P2 is the plane to the right of S. The function h = h(x) is called 
the scattering function of the instrument S. The instrument is called isoplanatic if 
h(x, y) = h(x — y). The image of an object f(x), located in the plane Pi is 

u(x) = h(x- y)f(y)dy, x G P 2 . 
J Pi 

The problem of the resolution ability of S in Rayleigh's formulation is stated as 
follows: 

Given two signals, fi = 5(x — a), $2 = 5(x + a), where S(x) is the delta-function 
and a is a vector, \a\ is small, will the function h ( x + a )+ h ( x ~ a ) }i ave one or t wo max- 
ima near x = 0? If two maxima can be detected, then, according to Rayleigh, the 
instrument S can discriminate between the two point sources the distance between 
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which is 2a. 

This criterion of the resolution ability is very limited: the class of objects to be 
discriminated consists of two point sources. Rayleigh's conclusion is : one cannot 
discriminate between two point sources if a ~ |, because the width of the function 
h(x) is of the order of A, where A is the wavelength. 

In [5]- [12] the resolution ability problem was examined in detail. In [5] a new 
definition of resolution ability of a linear instrument is given. This definition allows 
for a very wide class of the input signals, not only for two point sources (bright 
points), as in Rayleigh's definition. It takes into account the noise on the image 
plane, and it takes into account the properties of the procedure that reconstructs 
the input signal from its noisy image. The principal conclusion in [8] is : 

The resolution ability in Rayleigh's sense can be increased without a limit by 
apodization. 

In other words, the scattering function of a linear isoplanatic optical instrument 
can be made by apodization as close as one wishes to a delta-function on any fixed 
finite subdomain of the image plane. This conclusion is justified as follows. The 
scattering function is h(x) = j D j(y)e~ tx ' y dy, where x-y is the dot product of vectors, 
and D is a finite region, it is the output pupil of the optical instrument S, (see [3], 
Sec. 5.3), and j(y) is the field distribution on the output pupil D. If one puts a mask 
on the output pupil D, then the field j(y) is transformed into j(y)n(y), where n(y) 
is the refraction coefficient of the mask. Putting the mask on D results in the change 
of h: the new h = h n := f D j(y)n(y)e~ lxy dy. It is proved in [5], [6], [TJ, that, given a 
bounded domain A on the image plane, one can find a sequence of functions f m (y) 
such that the sequence S m (x) := f D f m (y)e~ iX ' y dy is a delta-sequence on A, that is 
linim^oo J A 5 m (x — s)(p(s)ds = 4>(x) for any smooth function cj) vanishing outside A. 
Moreover, the sequence f m (y) is given explicitly, analytically, in [6] and in [IT] , see 
also [HI p. 211]. Therefore one can choose n m (y) so that j(y)n m (y) = f m (y), and 
then the correspionding scattering functions h m (x) form a delta-sequence on the 
finite domain A, so that the width of h m (x) can be made as small as one wishes. 

By the Rayleigh criterion, this means that, for sufficiently large m the resolution 
ability (in Rayleigh's sense) can be increased without a limit. In this argument we 
did not take into account the noise in the image. It is proved in [5] that it is possible 
to increase the resolution ability (in the Rayleigh's sense) of a linear isoplanatic op- 
tical system without a limit by apodization even in the presence of noisy background 
in the image plane. The noise is assumed to be independent identically distributed 
with zero mean value and an arbitrary large but fixed finite variance. Our arguments 
are valid under the assumption that the classical Maxwell's equations are valid. We 
do not discuss quantum effects in this paper. 

The author thanks Prof. O.L. Weaver for discussions. 
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